Dictionary between scattering matrix and Keldysh formalisms for quantum transport 

driven by time-periodic fields. 
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We present the relation between the Floquet scattering matrix and the non-equilibrium Green's 
function formalisms to transport theory in noninteracting electronic systems in contact to reservoirs 
and driven by time-periodic fields. We present a translation formula that expresses the Floquet 
scattering matrix in terms of a Fourier transform of the retarded Green's function. We prove that 
such representation satisfies the fundamental identities of transport theory. We also present the 
"adiabatic" approximation to the dc-current in the language of the Keldysh formalism. 

PACS numbers: 72.10.-d,73.23.-b,73.63.-b 



I. INTRODUCTION 

In the last years, there has been an increasing theoreti- 
cal and experimental activity around quantum transport 
phenomena induced by time-dependent fields. Pumping 
phenomena in mesoscopic systems constitutes a very in- 
teresting case, where periodic out-of-phase potentials de- 
form the gates of semiconductor structures allowing for 
the generation of dc-currents even in the absence of a 
static bias, 1 . 2 . 3 . 4 . 5 . 6 . 7 ! 8 ! 9 ! 10 . 11 ! 12 ! 13 ! 14 

The scattering matrix approach and Keldysh non- 
equilibrium Green's function technique are the most pow- 
erful formalisms in the theory of quantum transport. 
Recently, the generalization of the Scattering (S) Ma- 
trix theory to time-periodic transport phenomena in an 
energy representation has been formulated 1 ^* 1 ^ 7 ", while 
an alternative treatment in a time representation has 
been proposed in Refs. Il8ll9t Keldysh formalism as 
a theoretical tool to investigate time-dependent trans- 
port phenomena in mesoscopic systems has been in- 



troduced some time 



,20.21 



and has been employed 



to study problems like ac-transport through quantum 
dots^S and superlattices 23 , Josephson junctions 2 ^*! and 
quantum pumps 2 ^ 7 .. Recently a practical formulation 
to treat problems with harmonically time-dependent po- 
tentials has been presented and used to investigate quan- 
tum transport in a mesoscopic ring threaded by a time- 
dependent nuji 2 ^ and systems with ac-potentials 29 ' 30 . 
Other formalisms to describe quantum transport in the 
presence of time-periodic fields are based in a modified 
transfer matrix approach 31 and in the Floquet represen- 
tation of the Hamiltonian with the introduction of non- 
hermitian dynamics for the wave function propagation, 
in order to represent dissipative effects gLM^ 

The Scattering Matrix formalism is basically a single- 
particle approach. Therefore, it cannot be directly ap- 
plied to systems described by Hamiltonians containing 
many-particle interactions. The theoretical framework in 



which Keldysh formalism is based is exactly the opposite 
one, namely the systematic treatment of many-particle 
interacting systems. This formalism is, however, also 
adequate to investigate transport phenomena through 
mesoscopic systems even in the case that many-body in- 
teractions do not play a relevant role. The reason it that 
the effect of the environment, in particular, the leads 
and reservoirs to which the mesoscopic system is con- 
nected, are suitably represented in terms of self-energies. 
In the description of quantum transport in systems of 
non-interacting electrons, the agreement between both 
formalisms is expected to be the rule. In the context 
of stationary transport, the equivalence between the two 
approaches was first pointed out by Fisher and Leei^ 

An important experimental situation corresponds to 
the case of slowly oscillating driving fields. The low fre- 
quency regime is sometimes loosely referred to as "adi- 
abatic" . This word stems from the Greek word "a- 
diabatos" , which means "not passable" . Traditionally, in 
theoretical physics, this term is employed when an iso- 
lated or closed quantum mechanical system is perturbed 
by a time-dependent Hamiltonian in such a way that the 
eigenstates do not mix as time evolves. This idea can- 
not be trivially exported to describe a quantum system 
coupled to an environment where the spectrum is contin- 
uum and concepts like energy levels are not well defined. 
In the framework of open quantum systems, the term 
"adiabatic" is sometimes understood as a synonymous 
of low frequency behavior while it is also sometimes em- 
ployed to define a description where the variable t in the 
time-dependent piece of the Hamiltonian is considered 
as a frozen parameter. An important number of works 
have been devoted to investigate quantum transport in 
pumps within the low-frequency regime. In particular, 
an "adiabatic" approximation to the Floquet scattering 
matrix has been introduced^S^iS'iS'i 7 " which, in prac- 
tice, allows for the evaluation of the contribution to the 
pumped dc-current that behaves linearly in the driving 
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frequency. 

The aim of this work is twofold. On one hand, we 
show that, for non-interacting quantum systems driven 
by time-periodic fields, in contact to static reservoirs with 
arbitrary densities of states, or with oscillating reservoirs 
described by smooth densities of states, it is possible to 
establish a transparent and complete dictionary between 
the Floquet scattering matrix approach of Refs. ll5llfl 17l 
and the Keldysh Green's function treatment of Ref. [2SJ- 
The second goal of this work is to formulate an "adia- 
batic" approximation, analogous to the one used in scat- 
tering matrix formalism, in the language of nonequilib- 
rium Green's functions. 

The work is organized as follows. In section II we sum- 
marize the description of quantum transport within the 
scattering matrix approach as well as the "adiabatic" ap- 
proximation to the Floquet scattering matrix. The de- 
scription of quantum transport for non-interacting elec- 
trons driven by time-periodic fields in the framework 
of Keldysh formalism is summarized in section III. We 
present several equivalent equations to evaluate the dc- 
component of the current flowing between the reservoirs 
and the central mesoscopic driven system. We consider 
cases where the pumping voltages are locally applied only 
at the central system, as well as cases where the voltages 
are applied at the reservoirs. A formula that allows for 
the translation between the two formalisms is presented 
in section IV. In section V we show that a fundamen- 
tal property like the unitarity of the Floquet S-matrix 
can be proved in terms of identities satisfied by the non- 
equilibrium Green's functions. In section VI we present 
the adiabatic approximation formulated in the language 
of non-equilibrium Green's functions. Finally, section 
VII is devoted to summary and conclusions. We have 
also included appendices A, B, C, D and E with relevant 
identities and properties used along the work. 



II. SCATTERING MATRIX FORMALISM 
A. General formalism 

The scattering matrix approach to quantum trans- 
port, usually referred to as the Landauer-Buttiker 
approachj2&2£ considers the propagation of carriers 
through a mesoscopic sample as a scattering process. The 
sample (scatterer) is assumed to be connected to several 
N r contacts (playing the role of electron reservoirs) via 
single channel leads. Then, all the information about 
transport properties of a mesoscopic sample is encoded in 
the scattering matrix S whose elements are quantum me- 
chanical amplitudes for electrons coming from some lead 
to be scattered into the same or any other lead. These 
amplitudes are normalized in such a way that their square 
define the corresponding particle fluxes (currents). 

In each lead there are two kinds of states, incom- 
ing to and outgoing off the scatterer. Correspondingly, 
we introduce two kinds of second-quantization operators 



a a (E) I a) a (E) and b a (E) /w a {E) which annihilate/create 
one particle in the lead a with energy E. Using these op- 
erators one can calculate a time-dependent current flow- 
ing into lead a as follows (we use a convention that the 
current directed from the scatterer towards the reservoir 
is positive) i^i 

I a (t) = f JJdEdE'e*^* 

(1) 

x{{bUE)b a (E'))-{aUE)a a (E')}}. 

Here (...) means quantum-statistical averaging over 
equilibrium states of all the reservoirs which are assumed 
to be unaffected by the coupling to the scatterer. 

1. Time-periodic local potentials and static reservoirs 

The particles incoming from some reservoir a are equi- 
librium particles. Therefore, if the reservoirs are station- 
ary we have: 

(aUE)ap(E')) = S a0 5(E - E')f a (E), (2) 

where f a {E) is the Fermi distribution function for elec- 
trons in reservoir a. In general, the reservoirs have dif- 
ferent chemical potentials \i a and temperatures T a . The 
oscillating potentials at reservoirs will be considered sep- 
arately. 

The operators b a for particles outgoing (i.e., scattered) 
into the lead a are related to the operators ap for in- 
coming particles (f3 = 1, . . . , N r ) through the scattering 
matrix^! In the present paper we will consider the scat- 
terer which is driven by external forces that are periodic 
in time with period t = 27r/f2 . Interacting with such a 
scatterer an electron can gain or loss some energy quanta 
nhflo, n = 0, ±1, . . . Therefore, in this case the elements 
SF,a/3(E n , E) of the scattering matrix (the Floquet scat- 
tering matrix Sf, see, e.g., Ref ill are photon-assisted 
amplitudes (times y/k n /k, k — y / 2mE/h 2 ) for an elec- 
tron with energy E entering the scatterer through lead 
(3 and to leave the scatterer with energy E n = E + uMIq 
through lead a. Then the desired relation between oper- 
ators b for outgoing particles and a for incoming particles 
reads as follows^ 

N r 

b a (E) =J2J2 S *M E , E n )a p {E n ). (3) 

13=1 n 

The sum over n runs over those n for which E n > 0. In 
general, some elements of the Floquet scattering matrix 
describe transitions between the bound states (E n < 0) 
and propagating (i.e., current-carrying) states or vice 
versa. Such processes do not contribute to Eq.JHJ). There- 
fore, in what follows, by the Floquet scattering matrix we 
will mean the submatrix corresponding to transitions be- 
tween propagating states only. In fact, if HQq *C E, the 
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sum in Eq.© runs over all the integers: — oo < n < oo. 
Note that if the scatterer is stationary, then only the term 
with n — remains non-vanishing. 

Current conservation implies that the Floquet scatter- 
ing matrix is a unitary matrixii£*i£ 



N r oo 

2_j 2_, S*p a {E n ,E)S l 3~ l {E n ,E ri 

(3=1 n=~oo 



(4a) 



N, 



^ ^ Sap(E, E n )S 1 p{E m , E n ) — S m o5 ai . (4b) 

0=1 n= — oo 

One can easily check that the unitary conditions guar- 
antee that the operators for outgoing particles obey the 
same anti-commutation relations as the operators for in- 
coming particles: 



[ai l (E),ap(E')]=5 a0 5(E-E'), 



[bi{E)ME')] = 5 a pS{E-E>), 



(5) 



However, in contrast to incoming particles, the scat- 
tered ones are non-equilibrium particles with the distri- 
bution function f { a out) (E)6(E - E')6 a p = (fit (£?)&+(£')) 
being different from the Fermi distribution function^ 



N r 



fL° at) (E) = J2 E \SFME,E n )\ 2 f a (E n ). (6) 

13=1 n= — oo 

Using Eqs.J2J and one can calculate the time- 
dependent current I a {t), Eq.©, flowing into the lead a. 
The dc component I a of this current reads as follows: 



/,. - / 'IE 

o 



{/("«>(£?) -/„(£?)}. (7) 



Substituting Eq.© into the above equation, using 
Eq. (|4bjl and making the shift E — ► E — riMl^ wc finally 
get: 



i« = f;jdEY: e i^^.^r 



/3— 1 n= — oo 



(8) 



x{/#)-/ a (a^ 



This expression emphasizes that for weak driving ampli- 
tudes and small frequencies, only electrons close to the 
Fermi level do contribute to the current. 

An alternative expression for the above current may 
be obtained by making the energy shift in /i {E) only 
and using Ea. (|4al) . We get a dc current in a more usual 
form (as a difference of forward and back photon-assisted 
transmission probabilities) : 



oo N r oo 

I a = ^JdEY l E { \SfMEu, E)\ 2 fp(E) 

j3=l n=-ca 

-\S F . Pa (E n ,E)\ 2 f a (E)}. 



(9) 



Given the above equations either JHJ or @ define a dc 
current flowing through the scatterer coupled to station- 
ary equilibrium reservoirs. 



2. Time- dependent voltages at the reservoirs 

If the reservoirs are subject to oscillating voltages with 
the same frequency: 



V a {t) = V a COs(Q t + <Pa), 



(10) 



then we proceed as follows lii^S Within the reservoir with 
uniform oscillating potential the electron wave function 
has the structure of the Floquet function type: 



-iEt/h-ih,- 1 f dt'eV a (t') 

y a (V a ,E,r) =e -°° ip E (r) 



iEt/ 



oo , \ 

V*M E Jn(^) 
n — — oo v 7 



(ii) 



We introduce operators a 1 ' (E) / a' (E) which cre- 
ate/annihilate one particle in the given above Floquet 
state. The corresponding distribution function is the 
Fermi distribution function: (a'£(E)a'g(E')) — S a/ 3S(E — 
E')f a (E), where the Floquet quasi-energy E is chosen 
to be equal to the energy of the stationary state with 
the same spatial part ipE (r). This follows from the fact 
that the uniform oscillating potential does not change 
the normalization of the wave function: J d 3 r\^!(E)\ 2 = 
/ d?r\ij}E\ 2 = 1- As a consequence, the only occupied 
states ^>(E) are those that correspond to the ones in the 
stationary reservoir with the same ipE- 

Then we construct operators a a (E) / a^E) corre- 
sponding to incoming particles in the leads with definite 
energy E as follows: 



a a (E) = 



n=—oo 



./„ (^)< iE-nh^). (1.2) 



The operators b for scattered particles are related to the 
operators a for incoming particles through the scattering 
matrix of the sample, Eq.©. Substituting the calculated 
a and b operators into Eq. © and averaging over the time- 
period we get a dc current in the presence of oscillating 
potentials V a (t) at reservoirs as follows: 

oo oo C N r 

i a = fJdE e E ME~nhn ) 

n=-oo L /3=1 



x E S Faf3 (E,E q )SF, a /3(E,E m ) 

m,g— — oo 



(13) 



x J, 
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Using Ea. i|4b[l and making the shift E — > E + nhilo in the 
last term one can rewrite the above equation in a more 
convenient form: 



JV. 



I a = iJdE £ E \fp(E-nMl )-f a {E)\ 

n=-oo/3=l > 



E Sp a0 (E,E q )SF,a/3(E,E m ) 



stationary scattering matrix and have to be calculated in- 
dependently, see Refll7l for simple examples. Notice that 
the current conservation introduces some constraints to 
the matrix A. Substituting Eq.lJTSJ into Eq.|@} and tak- 
ing into account that the stationary scattering matrix is 
unitary we get the following^ 



ih ( dS^dSo _ df^dSo 
~2 [ ~df~dE ~ ~dE~dt 



(17) 



f eV \ j ( eVy 

\ Mi a J Jn + m \ nn 



i(q-m)ipj3 



(14) 

To calculate the current flowing through the driven 
mesoscopic sample it is necessary to know the Floquet 
scattering matrix Sf- At strong driving this matrix has 
an infinite number of elements. Therefore, its calculation 
is, in practice, a non-trivial problem. This problem can 
be greatly simplified if the driving frequency is small. In 
this limit the Floquet scattering matrix Sf can be related 
to the stationary scattering matrix So having much less 
number of elements. 



B. Adiabatic approximation 

At low driving frequencies, Oo — > 0, one can expand 
the elements of the Floquet scattering matrix in powers 
of SIq- The lowest-order terms read as followstii 



S F (E n ,E) = S , n (E) + n f° g%fg) 

+hn A n (E) + o(n 2 ), 



S F {E,E n ) = So.-n(E) + n f° aSo -gZ lE) 

+nn A- n (E) + o(rt 2 ). 



(15a) 



(15b) 



Here Sq^ is the Fourier transform for the frozen scatter- 
ing matrix which is defined as follows: 



§o(E,t)= J2 ^ mnot SoAE)- 



(16) 



n— — oo 



Let the stationary scattering matrix Sq(E) depend on 
some parameters pi £ {P}, i = 1,2, . . . ,N P which are 
varied under an external drive. Since we assume that 
the latter is periodic, the parameters are periodic in time 
as well, pi(t + To) = Pi(t). The frozen scattering matrix 
So(E,t) is defined as the stationary scattering matrix 
So(E,{P}) with parameters being dependent on time: 
S (E,t) — S (E, {P(t)}). We emphasize that the frozen 
scattering matrix does not define the scattering proper- 
ties of a driven scatterer. Only the Floquet scattering 
matrix does. 

The matrix A(E,t), whose Fourier elements A n (E) 
that enters the expansion l|15fl . cannot be related to the 



The matrix A reflects a directional asymmetry of a dy- 
namical scattering process arising as a result of interfer- 
ence of photon- assisted scattering amplitudes. 40 

Equations i|15fl and (|17fl show that the expansion in 
powers of Qq is, in fact, an expansion in powers of 
hfto/SE, where SE is an energy scale characteristic for 
the stationary scattering matrix. The energy scale SE re- 
lates to the inverse time spent by an electron inside the 
scattering region (the dwell time). Therefore, alterna- 
tively one can say that the adiabatic expansion given in 
Eq. (|15|) , can be applied if the period of an external drive 
is large compared with the dwell time. This definition 
of an adiabatic regime is different from a usually used 
in quantum mechanics one which relates the excitation 
quantum M7o to the electron level spacing. 

We conclude, to find the Floquet scattering matrix 
with accuracy of order Oo it is necessary to calculate 
two matrices, So(E,t) and A(E,t), each of them having 
N r x N r elements. If this is done we can calculate the 
current up to terms of order Oo. For the simplicity we 
suppose that all the reservoirs have the same chemical 
potentials and temperatures, hence / Q (w) = /o( w ),Va. 
Then substituting Ea. ((r3|l into Ea. ((T4fl . expanding the 
difference of Fermi functions in powers of Oo, performing 
inverse Fourier transformation, and keeping at the end 
only linear in Oo and V a (t) terms we get the dc current 
as a sum of three contributions: 



I a = dE 



dfo(E) 
dE 



) / ^{i { r mp \E,t) + 



l£ ect \E,t)+lg nt) (E,t)}, 



(18a) 



8Sl(E,t) 



(18b) 



2 N r 

li rect \E,t) = ^J2{Vf,(t)-V a (t)}\SoME,t) 



13=1 



(18c) 



N, 



li mt \E,t) = f E Vp(t)l 2hn o Re[S* a0 (E,t)A aP (E,t)] 

13=1 



' 2 I dt dE dE dt 



(18d) 
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Here, in Ea. (|18djl . we have integrated over energy the 
term proportional to d 2 f /dE 2 and made it to be pro- 
portional to dfo/dE. 



The first contribution, 7, 



(pump) 



(E,t), is an adiabatic 



current pumped by a dynamical scattered This current 
is non zero if the time-reversal symmetry in the system is 
broken by the external drive, So(t) ^ So(—t). To this end 
at least two parameters of the scatterer have to be varied 
with a phase lag different from zero and n. The second 
contribution, Ia ect ^ (E, t), is due to the rectification of ac 
currents flowing under the influence of ac voltages by the 
varying conductance of a sampled! This contribution is 
non-zero if the potentials V a (t) are different. The third 
contribution, Ia nt \E,t), is due to interference between 
the ac currents generated by the dynamical scatterer and 
the ac currents produced by the external voltages^ This 
current can be non-zero even if all the potentials V a (t) 
are the same. Formally, this contribution can be viewed 
as due to external voltages acting as additional pumping 
parameters. So, if all the potentials are the same, V a (t) = 
V(t),Wa and only one parameter of a scatterer is varied 
then only the third contribution remains. 

Now, we turn to the Keldysh formalism to quantum 
transport for harmonically driven systems and then, we 
will establish a correspondence between the two for- 
malisms. 



III. QUANTUM TRANSPORT WITHIN 
KELDYSH FORMALISM 



A. Hamiltonian, Green's functions and Dyson 
equations 

The goal of Keldysh formalism is the evaluation of the 
Green's function of the system in the real time axis. The 
starting point is a Hamiltonian to describe the driven 
system in contact to the N r particle reservoirs through 
connecting leads: 



H{t) = H sys (t)+H c 



H res (t). 



(19) 



Although this formalism provides a systematic way to 
approximately treat many-body interactions, we focus in 
Hamiltonians corresponding to non-interacting electrons, 
which can be treated exactly. In several problems, it is 
convenient to describe the central driven system by a 
lattice model containing N sites. We adopt here that 
point of view. For such a system a generic Hamiltonian 
of spinless noninteracting electrons reads 



N 

H«"(t)= y j {H s l J{t)c\c v +H.c:, 



Ll' = l 



(20) 



We assume that the matrix elements contain static and 
time-dependent pieces of the form: iJ^, s (i) = £y/(<I>) + 
Viji (t, 5ij>), being $ a static magnetic flux and the terms 
depending on the driving fields being periodic in time, 



Vi,v{t,8 l>v ) = !X-oo e ~ lfcnot ^M'(fc) with amplitudes 
depending on the phases <5y/ and Vi t i>(0) = 0. For the 
moment, we do not write explicitly the dependence of the 
matrix elements of H sys on the magnetic flux and the 
phases. Instead, we simply write eij> and Vi^(t). We 
shall recover the more complete notation in section VI, 
where we shall analyze symmetry properties that depend 
on those parameters. In order to simplify the notation 
we also adopt energy units with h = 1 . 

The contacts between the system and the reservoirs 
are described by hopping terms of the form 



H cont = -J2w a (c{ m c ja + H.c. 



(21) 



where k a and j a are, respectively, coordinates of the 
reservoirs and the central system. We assume models 
of free electrons for the reservoirs, 



H ^ =J2e k J ka c ko 



(22) 



We also consider the possibility of applying ac-external 
voltages at the reservoirs, which are represented by 
potentials of the form l|l()jl. This introduces a time- 
dependent shift in the energies Ek a — > s a (t) — 
£fc„ + eV^cos^o^ + fa)- Alternatively, it is possi- 
ble to get rid of these time-dependent shifts by re- 
course to a gauge transformation Ck a — > Ck a (t) — 

Cfc Q e Jt o ' . Withm such a procedure, 

the hopping matrix element of 121|) becomes modified to 

—i f t dt\eV a cos(fioti+w a ) n j-i c ± i r 

w a — > w a e 'o . Both ways of tackhng 

the problem end in the same results for the evaluated 
mean values of observables. We choose the first one. 

The elementary theoretical tools of Keldysh formalism 
are the retarded, lesser, advanced, and bigger Green's 
functions. The first ones are defined as follows 

Gfy (t, t') = Q(t- t') \G> (t, if) - G< (t, t')} , (23a) 



Gf l ,(t,t')=i(cl(t')ci(t)), 



(23b) 



while the remaining functions are defined from the re- 
lations: G>,(t,t') = [GJ?,(f,t)]* and G£,,(t,f) = 
[G^(t',t)]*. These Green's functions are used in the 
calculation of mean values of observables, in particular, 
the currents flowing through the different pieces of the 
driven system. Their evaluation needs the solution of 
the equations governing their time-evolution, which are 
derived from a matricial Dyson equationa£. If we focus 
on spatial coordinates lying on the lattice of the central 
system it is convenient to define the matrices H sys (t) 
and G R,< (t, t'), which have matrix elements Hf y „(f) and 

Gilf" (t,t'), respectively. The corresponding matrices for 
the advanced and the bigger Green's functions are ob- 
tained from the relations G A {t,t') = [G R {t' ,i)]t and 



6 



G > (t,t') = [G < (t',<)] t . The Dyson equations for these 
Green's functions are: 



f* dh± R (t, ti)G*(ti, tf) = is(t - f), 
G R {t,t'){ -i^ -H s y s (t')\ 

f t ,dtxG R {t,ti)± R {t u t') = \5{t-t'), 



(24a) 



{ ? | _ H*vs(t) }6<(t, t>) - dht R (t, t 1 )G<(t 1 ,t') 
-f^ oo dt 1 £<(t,t 1 )G A (t 1 ,t>)=6, 
£<(t,f){-»jr -H s * s (t>)} - f^dhG^t, tx)t< (tut') 
-f* oo dt 1 G<(t,t 1 )£ A (t 1 ,f) = b > 

(24b) 

where d and d indicates that the operator acts to the 
right and to the left, respectively. The Dyson equations 
for the Green's functions with spatial coordinates outside 
the central region have a similar structure, considering 
the corresponding pieces of the Hamiltonian H instead of 
H sys (t) and £ R >>(t, t') = 0. In the absence of many-body 
interactions, the self-energies entering l)24al) and l|24b|) 
take into account only the effect of the "escape to the 
leads" . They are obtained from the Dyson equations for 
Green's functions with coordinates along the connections 
to the leads and integrating out the degrees of freedom 
related to the reservoirs ,22i2I Explicitly, they read: 



and 



E* (t, t') =Q(t- t') E> , (t, - E< , (t, 



(25) 



(26) 



with 



, — "0 "o 



■ik<p 



duj 
— < 
2tt 



-i U (t-t') x > ( w _ m o ) Pct ( w _ TO n ), (27) 



where the density of states p a (u>) — J2k ^i 1 ^ ~ £ k a )i cor ~ 
responds to H res , while A<(w) = f a {uj) and A>(w) = 
1 - f a (u), being f a (u) = l/( e M^l) + i) ; the Fermi 
function corresponding to the reservoir a, which is as- 
sumed to be at the temperature 1//3 Q . This function is 
only well defined for equilibrium problems and is thus 
related only to the stationary component of the Hamil- 



tonian H res iSi We also define the functions 

r Q (fc,w) = 

xr°(w-mO ) 



v v 

_ "0 "0 



r>(fc,w) 



y y 

/ , Jm(-TT-)J m +k(-^ L ) 
"0 "0 



m— — oo 



xAq (u> — mf2o)r°(w — mfio), 
\w a \ 2 p a {iu). (28) 



In the practical solution of the problem, the strategy 
followed in Ref. I28H29I was to work with convenient in- 
tegral representations of the Dyson equations l|24a|) and 
ll24Tt : 

G R (t,t') = G°(t-t') + 

J dt 1 e' tkn ° tl G R {t,t 1 )V(k)G°{t 1 -t') + 

dw 



k— — oo 

00 

E 



d tl / dt a ~ e -<*no*i e -*-(*i-*») 

27T 



xG H (t,ti)X)(fc ) ci;)G (ti-t'). 



(29a) 



G>{t,t')= I dh [ dt 2 G R (tM)^>{tiM) 



xG A (t 2 ,t'), 



(29b) 



where ^ fc ' denotes X^o- ^or ^ ne coordinates (j a ,k a ) 
along the contact, it is convenient to work with 



GfJt,t') = -w a fL„db 1 G?.(t ) t 1 )g*(t 1 ,t), 



G1 



(30a) 



«(*,f) = -Wa{/I 0o #lG£ J jMl)<£(*l,t') 

+ /!' oo dt 1 G< Jct (t,i 1 )^(t 1 ,t')}, 

(30b) 

where we have defined Gf;%(t,t') = Ek a G ^ a ^)' 
while g R (t,t') = Q(t -l')[g>(t,t') -l<(t,'t')} and 
g A (t,t') = [g R (t',t)]*. In equations (|2^ajl and (|29b1 . we 
have also defined: 



E(fc,w) 



duj' f(fc,u/) 



(31) 



where T(k,uj) has matrix elements Ti t ii(k,u)) = 
6ij a 5i> j a T a (k,uj). The retarded Green's function G°(t — 
t') corresponds to the equilibrium problem defined by 
the static piece of H sys (t) (with matrix elements £i,v) 
dressed by the static component of the self-energy. In 
other words, it is the solution of 



[wi-e-E(0,u;)]G ,o H = i, 



(32) 
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being 



dt'G (t~t')e t ^+ lQ+ ^ t - t '\ (33) 



Notice that in the equation for the lesser (bigger) 
Green's function, we have dropped a term that depends 
on the solution of the homogeneous equation, which is 
only relevant in the description of transient behavior. 

According to Refs. we introduce the following 

Fourier transform for the retarded Green's function: 



G R (t,Lu) 



dt'G R (t,t')e i ^ +i0+ ^ t - t '\ (34) 



Transforming (|29a|) according to l|34|) results in the fol- 
lowing set of linear equations 

G R (t,iu) = G°(u;) + 

e - tknot G R {t,uj + kno)V(k)G a {oj) + 

k— — oo 
00 1 

e ~ tknot G R (t,uj + kn )t(k,uj)G (uj). (35) 

k— — 00 

Since the above equation is periodic in t with period tq = 
f2o/27r, it is possible to expand its solution in Fourier 
series: 



G R (t,w) 



(36) 



k= — 00 



Note that the above procedure takes care of causality. 
In fact, the transformation 134fl . defined with respect 
to the difference of the two times is the natural exten- 
sion to the transformation l|33l) for retarded Green's func- 
tions defined in text books for stationary problems, which 
ensures correct analytical properties of the transformed 
function. In addition, notice that (|36|) is a Fourier series, 
not a Fourier transformation, which reflects the fact that 
Dyson equation is periodic in the "observational time" 
t and so does the corresponding solution. The retarded 
Green's function can be calculated from the solution of 
the linear set <|35[1 . A convenient method for the direct 
evaluation of the Fourier components (|36|l is the renor- 
malization method of Ref. l30t 



B. Current through the leads 

The mean value of observables related to one-body op- 
erators can be directly expressed in terms of the lesser 
(or bigger) Green's function. In particular, the charge 
current flowing through the lead from the central system 
towards the reservoir a, can be written (in units of e/H) 
as: 



J*(t) = 



k a 



4,0 



= 2w a Re[Gf aa (t,t) 



(37) 



Notice that the above current is related to the current 
Q through: I a (t) = eJ a (t)/h. Taking into account the 
Dyson equation describing the contact between the cen- 
tral system and the lead l|30b|) this current can be written 
as: 



J a (t) = -2|w; Q | 2 Re| / dt 1 

J — DC 

+Gf atja (t,t 1 )g£(t 1 ,tj\}. 



G 



(38) 



Making use of the definition (|23a(l , we can also express 
(|38|l as follows: 



G; 



,(Mi) 



J a (t) = -2|w Q | 2 Re| / dti 

G t j a (*> *0] ff« (*i > *) + G l j a &h)g£{h,t)} ,(39) 

which, for the case of stationary reservoirs, simplifies to: 



J a (t) = 2Im{ dh /+~ %L e Mt-H) T a a {uj)x 



The above expression has an appealing form, since 
it depends on Boltzmann-like factors. In fact, 
F° (u) f a (ui)G> a j a (t>tl) i s related to the probability for 
a state in the lead a to be occupied times the probability 
for its closest site at the central structure, j a , to be un- 
occupied, while the term (1 — / Q (w))r° (w)G^ *i) is 
related to the probability of the opposite process to take 
place. 

In what follows, we focus in the dc-component of the 
current, which is defined as: 



1 pr 

Ji c = -\ dtJ a (t). 
T o Jo 



(41) 



1. Review of the stationary case 

Let us first consider Vi^(t) — V a (t) = and we review 
the procedure introduced by Caroli et al4£ Since in the 
stationary regime the Green's functions entering eq. Ij38(l 
depend only on the difference of times, it is possible to 
perform the usual Fourier transform in the variable t — t'. 
The result is: 



./„ = -2 „■„ -'Re<{ / g 



+Gf a , 3a 



(42) 

Then, using the following identities and definitions: 

Rf \ r A/ M* f°° duJ ' Pct{u') 
9a( w ) =ifa(u)Pa(w), (43) 
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as well as the Dyson equation: 



G 



N r 

£j» = Y. G LiM*U") G UM> (44) 



0-1 



Eq. (|42() can be written as: 
"°° diu 



21m 



^Muj)r°(uj)\Gf a>j0 (uj)\ 2 } 



(45) 



Let us note that in the present case G R (u>) = G°(lo) with 
±(0,uj) = t°(uj), being 



a — 1 



(46) 



Being an equilibrium Green's function, G°(uj) satisfies 
the following property (see appendix A for a proof): 



Im 



Nr 

QjM] = E G W>)*)^,i>)*- (47) 

/3=1 

Using it in 1)45(1 and recalling the definition of the self- 
energy (|27l) . we get the well known expression for the 
current: 

Nr f°° d 

J « = E / ^r r ° H r °MI G L,>)l 2 //sH-/aH ■ 

(48) 



2. Time-periodic local potentials with stationary reservoirs. 

Let us now consider the possibility of time-dependent 
terms in the Hamiltonian of the central system, but sta- 
tionary reservoirs, i.e V a (t) = 0. 

Using (|29b(l as well as the Fourier representation 136|) 
in Eq. (|38[1 we get the following expression for the dc- 
component of the current through the lead a: 



J" 



N 



t = IZo ^{ 2Im [^^( ^)] r °(^)/ Q (^) + 



/3=1 ft=-oo J 

(49) 

There are two additional equivalent expressions to the 
above dc-current, which correspond to two different rep- 

. The first one 



resentations of the term Im 



is obtained from the condition of the continuity of the 
current, as shown in appendix B. An alternative proof is 
given in appendix C. Substituting i|B2|) in l|49f) results in 



the following representation: 

N r oo 



Jt c = E E / S>< 

\g ja ^{k,w)\ 2 Yl{w)Yl{u> + m a )f^) 

\g j0da (k,u>)\ 2 Tl(u;)r$(u; + kn o )f a (u;)}. (50) 



{ 



The second additional representation corresponds to 
substituting the identity (|D4|) derived in appendix D into 
Pty. The result is: 



N r 



dm 



J " C = E E / ^r°( w )r°( w + ^„) 

x\G ja<j0 (k,uj)\ 2 [f p (cu) - f a {u + fcfio)] . (51) 

Any of the three representations l|49|) , l|5U|l and 151|l are 
equally valid to calculate the dc-current flowing through 
the lead a. The concrete evaluation implies the solution 
of the retarded Green's function from the Dyson equation 
with t(k,u) = 0. 



3. Time- dependent voltages at the reservoirs. 

Let us finally consider the more general case, where, in 
addition to the pumping potentials at the central struc- 
ture, ac voltages are applied at the reservoirs. 

We start from the definition of the time-dependent 
current through lead a (|39[) and substitute there Eqs. 
(|2"T|) and |29bJ). Then, we use the Fourier representa- 
tion of the retarded Green's function (|36|l and take the 
dc-component. The result is: 

N T oo . 

/9=lAs,g,p=-oo °° 

[T a (p, u + kflo)Fp(q,u) - T<(p,uj + kQo)Tp(q,u)] 
x Sj a d„ (k-q+ P ,uj + qtl Q )g* a J0 (k, u;)}. (52) 
Equivalently, this expression can also be written as: 

Nr OO 



^ dc = E E wf^(f^) 

p — 1 — oo 



xW^)J m (^)Re{ 

xlj^(w + (n — /c — m)f2o) + ( n — k — m)fio) 

-fa(u) Sj a ,j (k -q + p,u+(n-k + q)Q ) 

xg* aj0 (k,u; + (n-k)n o )y (53) 

For reservoirs with a smooth density of states such that 
r° (lo — mfio) ~ r a (o;) and r^(w — mf2o) ~ ^(w), the 



2^ 
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expression for the dc-current further simplifies. In fact, 
in such a case l|52[) reduces to: 



oo f 

J d a c =E £ Re{/^r>)i>) 

(3—1 k,q,m— — oc 



e nvp x 



fp(co - mn )J m+q (^-)J m (^-)g ja ^ (:t (k -q,u + <?fi ) 



(fc, w) - e*™°f a {u + (k- m )n )J m+q {^t) 



xj 7 

(54) 

where we have used the first summation formula for prod- 
ucts of Bessel functions given in appendix E. Performing 
a shift to — ► cj — fcfio in the second term of (|54|l and mak- 
ing use of the identity (|D4fl . it is also possible to recast 
the above expression as: 



x 



qm— — oo ^ 

E E C"w§w> 

/3— 1 — oo 

^j'cJjb + <ltt Q )g* a j0 (k, u)fp{u) - mClo)) 

-ie tqv " J ™+?(^)^m(^) (g, w) 

~^i Q jc( _ 9 ; w + qQ )jf a (uj - mO ) |. (55) 

In summary, expressions (I52|) and (|53|) define two equiv- 
alent ways to calculate the dc current through the lead 
a in the case of ac voltages at reservoirs with arbitrary 
densities of states, while (15411 and l(55|) are two equiva- 
lent representations of such current for reservoirs with 
smooth densities of states. 

In any of these cases, the evaluation of J^ c implies the 
solution of the complete set ()35[) . 



IV. TRANSLATION BETWEEN THE TWO 
FORMALISMS 

We propose the following translation between the Flo- 
quet S-matrix and the Green's functions: 



(56) 



i^Jv^(uj + mfi )r^(w + nQ, )Qj a j l3 (m — n,uj + nfio), 

with E m = lu + mf^o- 

It is important to note that this translation exactly 
recovers all the representations for the dc current derived 
in the framework of Floquet scattering matrix theory in 
the case of stationary reservoirs. In fact, translating H5Uf) 
and (|51|) according to 1)57(1 . leads to equations © and 
©, respectively. In addition, in eq. (|50|) it is possible to 



identify the transfer matrix formulation of Refs. |3lJ and 



T a p(E,E k ) - |0 ia j,(M)rr2(w)C(w + An o ), (57) 
being a ^ 13. 

For stationary problems, we should consider m = n = 
and Gj a j (k,uj) — ► Gj ai j p (u) m jSt i m which case, we 
recover an expression like the one proposed in Ref . l35t 

S° a0 (E)=S a ,p- 

In the case of reservoirs with ac-voltages described by 
wide-band models with smooth densities of states, such 
that T a (oj ± fcfio) ~ r Q (cj), the translation (|57Jl also 
transforms Ea. (|55j) into Eq. ljT3|) . 



V. PROOF OF THE UNITARY PROPERTY OF 

S F . 

In the S-matrix formalism, current conservation im- 
plies that the scattering matrix is unitary, see Eqs.Q. 
In what follows we show that, for static reservoirs, this 
property can be proved by using identities satisfied by 
the Green's functions as well as the translation formula 
(|57|l . In fact, let us use Ea. (|57l) inside the summations of 
the left hand side of Eq. JIEJ|: 



N r oo 

SF, a p(E,E n )SF,-fp(E m ,E n ) — 

p—1 n— — oo 

^/r°(w + mfl )T° (u>) \iGj y ,j a (m,Lj) - 

N r oo 

^Lj 7 ( _m > W + m ^°) ~ X! X! ^J't J/3 ( m _ U + n ^°) 

/3=1 n— — oo 



xT^(w + nO o )0£,, i/3 (-n, w + uO )] • 



(59) 



From the identity (|D3|> , it can be shown that the second 
term of the right hand side of the above equation vanishes 
identically, thus recovering equation (|4b|) . 

Similarly, using Ea. (|B7jl inside the summations of the 
left hand side of Eq. jta}: 



N T oo 

X/ S Fta p(E n , E)S F , a -r(E n) E m ) = 

a=l n— — oo 

VA":0 + ^r^w + mf7 )r°(w)^^ J - )3 (m,w) 

Af». oo 
a=l n— — 'Do 



(60) 



and from the identity i|G6(l , equation l|4a|) is recovered. 
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The same procedure can be followed to prove the uni- 
tarity of the S-matrix in the case of reservoirs with oscil- 
lating voltages at the reservoirs, provided the density of 
states of the reservoirs is smooth. 



VI. ADIABATIC APPROXIMATION WITHIN 
KELDYSH FORMALISM. 

A. Definition 

The adiabatic point of view is inspired in a paramet- 
rical representation of the time dependent terms of the 
Hamiltonian. This means a description where the obser- 
vation time t is assumed to be frozen in the equations 
governing the dynamics of the system. In particular, in- 
stead of the Dyson equation (|35(l , in a frozen description 
we must consider the following equation 

OO 1 

k— — oo 
oo 1 

J2 e - iknot G f {t,u>)t(k,u;)G (u)) (61) 

k— — oc 

which is a stationary Dyson equation corresponding to 
the strength of the parameters V(t) and V a (t) at the 
observation time t. As the potentials are periodic in t, 
the frozen Green's function can be expanded in a Fourier 
series: 

oo 

G f (t,u>) = e- iknot g f (k,oj), (62) 

k— — oo 

and through the translation H58|) it is possible to define 
the Fourier coefficients §o(k,E) = So,k(E) for the ele- 
ments of the frozen S-matrix as 

So, a p(k,E) = 8 a ,p8k,o ~ i W^MT^M Gj a , jfj (k , u) . 

. , . (63) 

In the Floquet S-matrix formalism, the adiabatic ap- 
proximation is given by Eas. p5|l . In analogy to Ea. (|15a( l 
we propose the following ansatz for the oc Qq approxima- 
tion to the Green's function: 

g R (k, U ) ~ gf(k, U ) + ^£ ggfe^) + n a(k, u ), (64) 

2 ouj 

or, equivalcntly, 

G(t, w) ~ & (t, lo) + 4 ^g^^ + n oa(t, u>), (65) 

where Q*(k,uS) (or &(t,ui)) is the frozen Green's func- 
tion, which obeys the equilibrium Dyson equation (|61|l . 



The substitution of the ansatz Eg. (|64fl with the trans- 
lation formula (|57|l into Ea. l|15a[) leads to the following 
relation between A a p and a, 



A a0 (k,E) = -^ro( w )r»( w )a ja j,(k lW ) 

i 9r» i sr°M 



. kflo 



r°H du r°(w) duo 

A a0 (t,E) = -i^( w )r°( W )a jaJ? (t, W ) + 



i dr° a (u 



at 



(66) 



InRef.Ha 

some important properties of the matrix A 
have been proved on the basis of the unitary property 
of Sf and the fact that for a Hamiltonian of spinless 
fermions that depends on a magnetic flux (f> and on time- 
dependent potentials of the form 

Vi, v (t) = 5 ttV hi W? + Vi c os{ut + 5j)], (67) 



the stationary matrix transforms under t — > — t as: 

SoMk, E, $, 5) = SoM-k, E, $, -S), (68) 
as well as 

So^Qe, E, S) = SoM k > E, $, 5). (69) 

Analogously, the Hamiltonian is invariant under the 
simultaneous change of t — > — t and 8j — > —8j. Therefore, 



= glA-k^-s^w). (70) 

In addition, in the presence of a magnetic flux <E>, the 
static terms of the Hamiltonian of the system satisfy 
£M'($) = = fo',i(*)]*> which implies: 



GUt,*,w) = G{,,(t,-$,i 



(71) 



In other words, the frozen Green's function ^ (k,ui) 
has the same symmetry properties of So,a/3 There- 
fore, from Eq. I|66ll we see that dj a j l3 has the same sym- 
metry properties as A a p. 

In order to calculate a(k,u>) explicitly, we have to con- 
sider the Dyson equation for the Fourier coefficients of 
the retarded Green's function, which can be obtained by 
expanding l|35|) in Fourier series: 

g{k,w) = G°{u)5 k>Q + 



OO 

£ 

k' — — oo 



g(k + k',LU + k'n Q )V(k')G Q (cu) + 
J2 g(k + k',u} + k'n )t(k',cj)G ((j), (72) 



k' — — oo 



11 



then substitute (|64|l and keep terms up to the first order 
in fl . The solution of the ensuing linear set allows for the 
evaluation of a(k, oj) as a function of the stationary Green 
functions G°(u>), the frozen Green functions Q^(k,u)) and 
the derivatives dQ s (k,oj) / 'dui and dG°(uj)/duj. Alterna- 
tively, in order to get a(t,Lu), we have to substitute the 
ansatz l|65|) in the Dyson equation Ij35(l , keep terms up to 
ex fio and solve the resulting linear set, which gives a(t, uj) 
as a function of &(t,u), d 2 G? (t,u)/&tdu, G°H, and 
dG°{Lo)/duo. 



B. Calculation of the dc current 

1. Time-periodic local potentials with stationary reservoirs 

In order to calculate the adiabatic approximation to 
the dc-current through the leads in the case of station- 
ary reservoirs it is convenient to start from Eq. Q51JI. 
Expanding that expression in powers of f2 and keeping 
up to the linear term, the dc-current reads: 



jdc 



r Jo 



2tt 



(ftM - f a (uj) (j^ as) (t,co) + Ji mt) (t^))] , (73) 



being 



j^)((, w ) = r»Hr»Hx 



0=1 



dt 



+ c.c. 



, (74a) 



J^Ht,u) = ]T |G^(i, W )| 2 r>)I>), (74b) 
0=i 



N,. 



isr°(u,)_ 



0=1 



2 duo 



dt 



+ c.c. 



-<gU(*.«) 



dt 



i UJ )a ja j l3 (t : u}) + c.c.) 



(74c) 



For the case of reservoirs with identical chemical poten- 
tials and temperatures, such that f a (u>) = fo(cu), Va, 
only the pumping term (|74aj) contributes. This term is 
x Qo an d can be shown to be equivalent to Eq. l|18b|) 
through the translation formula l|57l) . The second term 
is the usual stationary contribution l|48|) . The "interfer- 
ence" term (x Clo([ip — Ha)) is a small contribution in the 
limit of small static bias Hp — /i Q . However, it may give 



rise to interesting behavior in the presence of magnetic 
field. In particular, in a two terminal set-up it is an odd 
function of a magnetic field ^i, in striking contrast with 
a stationary conductance which is an even function of a 
magnetic field. 



2. Time- dependent voltages at the reservoirs. 

In order to derive the oc fio contribution to the dc cur- 
rent, we substitute the adiabatic approximation to the 
Green's function Ij64(l in (|52J) and we expand the remain- 
ing terms up to 0(Hq). The latter step results in the 
following expansion: 

T a (p,uj + kfto)T^(q,u) - T<(p,uj + k£lo)Tp(q,u) 

j ( SY^\T (1Yel\t (-YB-\i ( eV < 3 \ 



m.n — oo 



xe ^+^){T a (uj)Tp(uj)[f p (uj) - f a (w) 



n 



(k - n)gi{w) - mg 2 {uj) 



g 3 (uj)(k - nf 



+m 2 g4,(u>) + m(n — k)g5(u>) | 



(75) 



being 



ffi(<*>) 



92 M 



03 M 



04 M 



95 (w) 



}. 



duj 



-T a (uj) 



-r a (w) 



dtu 2 
d 2 f a (u) 
duj 2 



fp{u) - f a (oj) 



} 



r a (w) ra 2 r^( w ) 



<9w 2 
d 2 fp{LQ) 
duj 2 

dT a (uj) dTp(uj) 



duo 



8lo 



//3(w) - /a.(w) 

_}■ 

//?(w) - / (w) 



(76) 



When substituting in (|52|l we use the properties of the 
Bessel function enunciated in appendix E and keep only 
terms proportional to f2o and V a . The resulting expres- 
sion shows a rather compact form in the case of reservoirs 
with smooth densities of states and the same chemical po- 
tentials and temperature, such that / Q (w) = /o(cj),Vq!, 
in which case gs(w) = 0, while gi(u>) — —g2{u) and 
g$(uj) = — <?4(w). Terms cx d 2 f (oj) / du; 2 can be reduced 
to terms oc d fn{uj) / du by integrating by parts. The final 
result can be written by collecting the different terms in 
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three kinds of contributions as in Eq. i|18|) : 

~tJ q M J_ 00 2n du> 
x j(P" m f) (t, u) + J { r ct) (t, u) + J ( a mt) (t, wj\ , (77) 

with Ji pump) (t,u) given by d^) and 



xT a (uj)Tp(uj)[v a (t) -V (t) 



(78a) 



v it){\ 

+2fi 



3aJ(3 ( 



dui dt 
GLj^^d^^r + c-c]}. (78b) 



Through the translation formula Ij57(l . we can identify 
the t hree terms (|7l"al) . (T7Hal) and i(78b"|l with l(TFE)l . (|T%c|) 
and i|18d|) . obtained with the S-matrix formalism. In 
the derivation of the expression for the interference term 
(I78b|l we have made use of (|D4(I expressed in the adi- 
abatic approximation (|64|l and an equivalent equation 
satisfied by the frozen Green function: 

CiJ fc > w )-C^( fc > w )* = 

E SL, ( k + k '> ^»^L , j0 ( k '> ">)* > (79) 

k'/3 

which leads to the following relation: 



It t r >)sLt(t + ^ L '^ 

/3=1 fe'=-oo 
N r oo 

= E E r >){sL ( k + k '^) a ^A k '^r 

13=1 k' = -oo 

+%.a,(*+'=>)sL»(< ! >>* 



'Jo J/3 



9w 



(80) 



for the derivation of several useful equations to calcu- 
late the current (in particular, the dc-component of the 
current) flowing between the reservoirs and the meso- 
scopic system. We have considered two different situ- 
ations: (i) Driving induced by voltages applied at the 
central structure, (ii) Driving induced by voltages ap- 
plied at the reservoirs. In both cases we have considered 
reservoirs with a general density of states. We have also 
proposed an expression that enables the translation be- 
tween Green's function formalism and the Floquet Scat- 
tering matrix formalism of Refs. In the case 
(i), we have shown that this formula is able to trans- 
late exactly the expressions for the dc current flowing 
through the leads obtained in the two formalisms. Fur- 
thermore, we have shown that it is enough to derive the 
unitary property of the scattering matrix by recourse to 
properties of the Green's functions. In the situation (ii) 
we were also able to translate expressions for the current 
and to demonstrate the unitary property of Sf if we as- 
sume within the Green's function formalism models of 
reservoirs with a smooth density of states. 

We have also formulated the so called adiabatic ap- 
proximation to the dc-current in the framework of the 
non-equilibrium Green's function formalism. We have 
used it to derive the different contributions to the dc- 
current linear in the driving frequency in the two situ- 
ations (i) and (ii) described above. Making use of the 
translation formula of section IV, it is possible to com- 
pare these expressions with the ones previously derived 
in the framework of Scattering matrix theorjii^iSiii. The 
equivalence is complete in the cases of stationary reser- 
voirs as well as in the case of oscillating reservoirs with 
a smooth density of states. 
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VII. SUMMARY AND CONCLUSIONS 

Starting from the formulation of Kcldysh non- 
equilibrium approach to systems in the presence of time- 
periodic fields of Ref. l29i we have shown several iden- 
tities satisfied the Green's function. This has allowed us 



APPENDIX A: 

In the stationary system with Vij/(t) = V a (t) = 0, the 
Dyson equation for the retarded Green function is simply 
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Eq. (|32() . From there, we can write 



G°(w) -G°(w) t = 
= {iw-e-E (w)} _1 - llw-^-SV) 1 } -1 

= {iw-e-E (w)} _1 

xjiw - e f - S^^tUiw - ft - E°( w )t}~ 

- (lu-i-tP{u)\ jiu;-£-E (a;)} 

xjiw-^-S^o;)^ -1 
= G°(lu){± (lu)-± (lu)^G°(lu)\ (Al) 

where we have used the fact that i = e' . 

APPENDIX B: 



Let us start from the definition of the dc current flow- 
ing through the lead a, Ea. (|49(l . The condition of the 
continuity of the current implies: 



E J t c = 



N, 



2 E / Sj~j.(o>") /«(«>) +E E 



a=l fc= — c 



da; 
2^ 



//3HI^^(fc^)l 2 r^Mr°( w + fcno). (bi) 



Since the above equation must hold for arbitrary chem- 
ical potentials and temperatures of the reservoirs, the 
terms multiplying the Fermi functions have to van- 

ish for each a, which means: 



-21m 



N r oo 

<7 iajQ (0,o;)] =E E |e wjQ (fc,c)| 2 rO( W +fctt ). 

0=1 k=-oo 

(B2) 



spacial indices of that equation. The result is: 
{i^ + (u/ + *0+) - H«"{t)}& R {t,u/) 

<a 1 e ib/ ^- t ^t R (t,t 1 )d(ti,u') = i, (ci) 

> 

G R (t, W )t{ + (a, - z0+) - H*«"(t)} 

dt 1 e- <t '( t - tl >(3' K (ti,w) t f: H (t,ti) t = i.(C2) 

Note that, unlike the stationary case, these equations are 
not simplified to a linear set of equations and the inverse 
of the Green's function must be represented in terms of 
integro-differential operators. By multiplying (|C1|I from 
the left by G R {t, w)t and g2) from the right by G R (t, lj') 
and then subtracting the two resulting equations we get 

G R (t, w )t - G R (t, J) = [G R (t, u>)^G R (t, to') 
(u'-w)G R (t,u^G R (t,u;')- 

dt 1 e luJ '^ tl} G R (t, u)t£ R (t, h)G R {tx, w') + 



— oo 

t 



dt ie -M*-*0 wjtiJ^i, * 1 )t^r«(t 1 to 1 ). (C3) 



If we now perform the expansions in Fourier series for 
G R {t,uj) and G R {t,uj')\ we obtain: 

0(-fc,w)t-0(fc,w') = 

OC 

(uZ-w + jfefio) E S(k',oj)^g(k + k',u;') - 



APPENDIX C: 

We present the proof of an important identity for the 
retarded Green's function. This identity has been previ- 
ously proved for the case of static reservoirs in Ref. |13 
within the framework of a different formalism. 

We first perform a Fourier transform in t — t' in the 
first equation of the set (|24a|l and take the adjoint in the 



E G f (k', cj){£{k", uj' + ik + k'- fc")Oo) - 

k' 7 k" — — oo 

E(k", uj + fe'fio)* }<3{k + k' - k", u/). (C4) 
For u>' = lo — fcfio, the above equation reduces to 

g(-k,ojy -g(k,u}-kn ) = 

E £(fc>)t{Et(fc'> + (fc' -k")Q )- 

k',k" 

±(k",uj + k'n )}g(k + k' - k",oj- fcfio). (C5) 

In the case of stationary reservoirs, the above equation 
further reduces to 

g{-k,Lo)^ -g(k,tu-kn ) = 

oo 

E ^(fc',w) t [E (a) + fc'^ ) t - 

k' — — oc 

E (w + fc'O )lg(fc + fc',w-fcO ), (C6) 
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which, for k = 0, reads: 

oc 
k— — oo 

£°(Lj + knoy}G(k,u). (C7) 



Note that the identity 1|B2|) derived in the previous ap- 
pendix is a particular case of this equation. 

The above identities are also valid in the case of oscil- 
lating reservoirs with smooth densities of states, as can 
be verified by using the fact that r° (lu + kilo) ~ T° (cj) 
in (|31|l as well as the first of the summation formulas of 
products of Bessel functions of appendix E. 



APPENDIX D: 

In this appendix we prove another important identity 
satisfied by the Green's function. We start from the defi- 
nitions l|23a|) and l|23b|) of the different Green's functions 
in Keldysh formalism and use Dyson equations for the 
lesser and bigger components (|29b|l . as well as the Fourier 
representation of the retarded Green's function. We get: 



N r oo roo , 



X 



a=lk 1 k 2 k 3 -oo J ~°° 

e _i [£l ,( t -i')+n (fc 1 t-fc3f')]g( fcl) u + k 2 n )g{k 3 , w)t. (Dl 



Transforming the above Green's function according to 
(|34fl and (|36[) allows us to write the following spectral 
representation: 



cm- E/^ 

k>k" = -oo J -°° Z7T 



Qjk + fc', J_ + k"Q )t(k", Lo')g(k>, tjQt 
X cj - (u)' + k'Cl ) + i0+ [ ' 



For the case of stationary reservoirs, it is easy to prove 
from the above expression the following identity: 



g(k,uj)-g{-k,uj + kn ^ = -% J2 G{k + k',uj-k'n ) 

k' — — oo 

xt°(w -k'Q )Q(k',0J -fe'Oo) 1 , (D3) 
which for k = reduces to: 



g(o,uj)-g{o,Lj^ = -i J2 G(k',u-k'n ) 

k' — — oo 

xt {aj-k'n Q )g(k',oj-k'n Q )l (D4) 



The above identities are also valid in the case of reservoirs 
with oscillating voltages provided that they are described 
by a wide-band model with a smooth density of states 
such that r°(w — mfl ) ~ r°(w). This can be easily 
proved by using the first of the summation formulae of 
appendix E. 

APPENDIX E: 



Products of Bessel functions satisfy the following sum- 
mation formulas: 



Jn+p{X)J n {X) — 5 Pt o, 

i— — oo 
oo 

^ J n +p(X)J n (X)n = X(6 Ptl + S Pt -i) 



J2 Jn+ P (X)J n {X)r 



X 2 



-$P,o ^ X 



n— — oo 



1 1 X 2 
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